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Abstract 

We prove that the stress tensor, r afe , of a molecular system with 
arbitrary, short-range interactions can be point-wisely expressed as 
the functional derivative of the partition function with respect to the 
local deformation tensor. In this approach the set of components of T ab 
has a simple interpretation as the set of Lagrangian multipliers which 
need to be introduced to enforce the conditions relating point particle 
displacements to the body local deformation tensor. The question of 
the non-uniqueness of the formula for r ab is discussed. 
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1 Introduction 



It has been shown in pQ that the invariance of the partition function of a molec- 
ular system under the change of variables induced by the (measure preserving) 
canonical transformation 

[{r }, {p }] - [{r '}, {p'}] , (1) 
C = r? + e a (f l ), i = l,2,...,N, (2) 
_ dr? (r ) „ _ r fe cV(r ) n 6 

^ r^"^ + cV r=rJ P ' W 

where e (r ) represents the infinitesimal displacement of the body elements H at the 
point r, implies for a system with short-range interactions the local equilibrium 
condition [2] 

VV fc (r)+^ xt (r)=0 (4) 

with T ab (f) the stress tensor at the point r and ^2ct(r ) an external force (density). 
The local expression of T ab (f) in terms of the degrees of freedom (dof's) of the 
elementary constituents of the body was derived by explicitly computing the func- 
tional derivative of the partition function with respect to the particle displacement 
and comparing with eq. ([4]). The generality of the approach guarantees its valid- 
ity in any statistical ensemble, for whichever type of (short-range) potential and 
boundary conditions, and in a classical as well in a quantum mechanical setting 
(see p~] for details). 

The purpose of this brief note is twofold. First of all, after recalling how the 
equilibrium condition ([!]) can be deduced from purely thermodynamic considera- 
tions, we prove that to the set of stress tensor components can be given an elegant 
interpretation as the set of Lagrange multipliers that are needed to enforce the re- 
lation between the displacement vector, e (r ), and the deformation tensor, r] ab (f), 
given by the formula [2] 

^(r) = i[V a e b (r) + V fe e a (f)]. (5) 

Secondly we rederive the local expression of the stress tensor for a molecular system 
(endowed with short-range interactions) by moving from the "passive" interpreta- 
tion of the transformation ([2]) (followed in pQ) where eqs. ([2]) and ([3]) are seen as 
a mere change of variables, to an "active" one where we imagine that e(r ) is the 

1 A sum over repeated spatial indices (a,b,... = 1, ...,£)) is understood. We also 
frequently use the shorthand notation df(r )/dr a — V a /(r ). 

2 Boundary conditions matching those imposed to the body must be obeyed by the 
function e(r). 
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actual infinitesimal displacement of the body elementary constituents at the point 
r. 

As a fall-out of the approach we present in this paper the question of the 
uniqueness of the stress tensor [6] can be neatly addressed with the conclusion 
that no ambiguity affects the formula (|23p below. We will show, in fact, that there 
is no freedom to add to this expression any arbitrary divergenceless, symmetric 
rank- two tensor, as a priori geometrically allowed by the structure of eq. 



2 Mechanics and thermodynamics 

Let us start with a discussion of the physics of local body deformations. Assuming 
that the system is at mechanical equilibrium at fixed temperature, the principle 
of virtual works [3] ensures us that the work done by the body under an infinites- 
imal local deformation is zero. Furthermore, if the deformation transformation is 
reversible, the variation of the free energy will be equal to minus the work done by 
the body [3] ■ In these circumstances one can write the variation of the free energy 
of the system under an infinitesimal (reversible) deformation, as a function of the 
particle displacement vector and deformation tensor, in the form 

dA( V ,e) = -6 rev L( V ,e) = f d 3 rr ab (f) V ab (f) - f d 3 r ^ xt (r)e a (r ) , (6) 

Jv Jv 

where the first term in the last equality corresponds to the work done by the body 
deformation and the second to the work done by the external force (if there is 
one). As recalled above, the sum of the two contribution vanishes provided the 
displacement vector and the deformation tensor are related as in eq. ([5|). 

The way to see what the condition of thermodynamic equilibrium implies 
for this constrained system, is to introduce Lagrange multipliers, X ab , to enforce 
eqs. ([5]), and define the "unconstrained" variation of A 



<IA. 

i ab 



fa,e;A) = [ d 3 rr ab (f) V ab (f)- I d 3 r ^ xt (f )e a (f ) + 
JV JV 



<f'r A""( r) ~rj ab (r) - ^[V a e b (f ) + V b e a (r )]] . (7) 



Imposing the vanishing of cL4 uncon (?7, e; A) immediately yields the relations 

T ab {r)-\ ab {r) =0, (8) 
^ xt (r) + V b X ab (f) = 0. (9) 

Eqs. ([5]) provide the announced interpretation of the set of stress tensor compo- 
nents as the set of Lagrange multipliers which enforce the constraints ([5]). Elimi- 
nating X ab between eqs. (|8|) and Q gives back the body equilibrium condition 
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3 Statistical Mechanics 



We now want to get an explicit expression of the stress tensor in terms of the 
elementary dof's of the system. To this end we need to display the functional 
dependence of the free energy on rj ab and e a . Upon comparing with the form of 
eq. ([6]), one can derive the desired formula for T ab (f). 

The procedure outlined above can be straightforwardly implemented in Statis- 
tical Mechanics. Let us consider, in fact, a system interacting through the short- 
range potential W[{g}] and let W ex t [{<?}] be a generic external potential. Working, 
for concreteness, in the canonical ensemble (but the argument that follows would 
similarly go through in the micro- canonical ensemble pQ), one has for the free 
energy the formulae 

A = -^logZ° c , (10) 

r N N 

Z c = J IK dD P) J v II(^) 6X P ( " /^cxtlM, {P}]) , (11) 

n° cxt [{q}, M] = n°[{ q }, M] +u cxt [{ q }} , (12) 

n°[{q}dP}]=^^+U[ti}}- (13) 

i=i 

In eq. the symbol Y\ N (d D p)U N (d D q) is a short-hand notation for the D- 
dimensional integration measure over the system phase space and V is the volume 
of the box in which the system is contained. The superscript "°" in the previous 
equations is to recall that they refer to the undeformed body in equilibrium. 

To find the functional dependence of the free energy upon rf b and e a , we 
have to provide the expression of the Hamiltonian of a system subjected to a local 
deformation (of the type indicated in eq. (|2]) ) . To this end we first notice that under 
the infinitesimal displacement e{q ) the line element squared changes according to 
the formula [2] 

dq a dq a -> dq a dq a + 2r] ab (q)dq a dq b , (14) 

with rj ab (q) related to e a (q) as in eq. ©. Consequently also the kinetic energy 
of the system will get modified because of the addition of the extra contribution 
coming from the second term in eq. (I14p . In fact, from eq. (114j) one formally gets 
for the modulus square of the velocity 

dq a dq a dq a dq a b dq a dq b 

dt dt ^ dt dt +ZV [q) dt dt • [ b) 
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The Hamiltonian of the deformed system will thus read 



H ext [{q},{phv,e] =H° ext [{q},{p}} + 

N a b i N N 

-E^)^r-o E ^(%)44[M]-E ea ®)^ex t [M], (16) 

1=1 J7"t=l l =l 



where we have introduced the definitions 

= g» - # , gij = Jtf . (18) 



The first term in the second line in the r.h.s. of eq. fjl6[> comes directly from eq. (|15p 
after passing from velocities to canonical momenta. The second and third term 
arise as a consequence of the particle displacement q a — > qf + e a (<fj) (eq. (|2|)) in 
^/[{g}] and U ex t [{<?}]> respectively. While the structure of the third term is obvious, 
the form of the second needs some explanation. First of all, we observe that we can 
always consider a translationally invariant interaction potential as a function of the 
set of the two-particle distances {qij}. Secondly, since U is assumed to be short- 
range on the macroscopic scale over which e can appreciably vary, we conclude that 
one can get non-vanishing contributions to the r.h.s. of eq. (|16p only from terms 
where all particle distances are very small (smaller than some typical microscopic 
length). In computing the variation of qij under a particle displacement, we are 
then justified in writing 

q? ~ q a j - q? + e a {q t ) - q] - e a {q 3 ) = V b e a (q t )(q b -<$ + •••, (19) 

9» - Hi + Uv a e b m + V^mf-^A. . . = qy + r) a \$)^A . . . , (20) 

qij Qij 

where dots represent terms of higher order in the differences qfj, that we neglect. 
When eq. ([20]) is introduced in W [{(/}], the second term immediately emerges by 
expanding in the small quantity r] ab . 

We stress that, as expected, in (fT6j) the body deformation is completely de- 
scribed by the tensor rj ab , while the displacement e a is directly coupled only to the 
external force. 

Inserting the Hamiltonian (|16p in the formulae for the partition function and 
free energy, one obtains the sought for functional dependence on rj ab and e a . To 
first order in 7] ab and e a one thus gets for the free energy variation 

N „ N 

-PKxtlivUp}} 



dA( V ,e) = ±J U(d D p) J v U( dD <l)e 
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r N n a n b 1 N N 

■ -E^(«)^ - o E ^(«)«&4Hff}] - E ea (*)^cx t [{9}] • (21) 

For an easy comparison to the equations of sect. [2] it is convenient to introduce in 
each term of the sum over the index i, the identity f v d s r 5(f — qi) = 1. Having 
done this, eq. (|2ip can be cast in the form 

dA(v, e) = 

r N 9- b i N 

= j rf 3 r/(r)(-^5(f-^-- E 8(r-®<$,J%[{q}]) + 
V i=l "' i^»=i 

r N 

- / d\e a {r )(£ 6<? ~ %)^ ext [M]> , (22) 

where (. . .) means ensemble average. Comparing with eq. ©, the expression of 
T ab (f) in terms of the elementary dof's of the system is readily identified as the 
tensor that multiplies the deformation tensor in the formula for the work done by 
the system under a local deformation. One finds in this way 

N a b i N 

r a V) = -(E^-^(^ + 2 £ (23) 

i=l j&i)=l 

This formula is in agreement with [1] and, once integrated over volume, with the 
expression that can be found in the classical papers of ref. [5]. 

In closing we observe that the approach we have developed allows to answer 
the old question of whether a stress tensor obeying eq. (jlj) is unique or not [6]. 
The question arises, as from a purely geometrical point of view, one could imagine 
to add to whichever expression of r ab an arbitrary diver genceless, symmetric rank- 
two tensor, still fulfilling eq. dH). However if, as we advocate in this paper, T ab 
is identified as the tensor which multiplies the deformation tensor, r] ab , in the 
formula which expresses the work done by the system under an infinitesimal local 
deformation (eq. ©), such a freedom does not exist anymore. 



4 Conclusions 

In this note we have given a consistent derivation of the microscopic expression of 
the stress tensor of a body, that complies with the principles of Thermodynamics 
and takes properly into account the geometrical constraints existing between the 
particle displacement vector and the body deformation tensor. In agreement with 
what is known to happen in the case of a continuum system, we find that the 
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possibility of defining a (local) stress tensor rests on the assumption that the 
interaction potential between the body elementary constituents (for the rest fully 
arbitrary) is "short-range". The discussion we give is totally general and holds in 
any ensemble and whichever the boundary conditions imposed to the system might 
be. Remarkably the whole procedure goes through also in a quantum-mechanical 
setting [T]. 

A consequence of the line of reasoning we have presented above is that appar- 
ently there is no room for any ambiguity in the expression of the stress tensor we 
derive, if identified with the tensor which multiplies the deformation ten- 

sor in the formula for the work done by the system under an infinitesimal local 
deformation. 
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